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Introduction. This work is motivated by a recent breathers existence 
proof in the one dimensional FPU system, given by the equations: 

X n = V'(x n+ i - X„) - V'(x n - X n -l) , tl € Z , (1) 

where V is a smooth interaction potential satisfying V'(0) = and V"(0) > 
0. Using a center manifold technique 2 , one can prove the existence of small 
amplitude breathers (SAB) with frequencies uib slightly above the phonon 
band if B = \V"{Q)V^{0) - (U( 3) (0)) 2 > 0, and their non-existence for 
B < 0. Our aim is to test numerically the range of validity of this theoretical 
result and to explore new phenomena. For this purpose we shall fix V(u) = 
u 2 /2 + au z + \u i , which yields B = 3(1 - 12a 2 ). 

We work with the difference variables u n — x n — x n —i more suitable for 
the use of our numerical method. We also use periodic boundary conditions 
Un+2p(t) — u n (t) so that the maximum frequency of the linear phonons is 
exactly 2 as in the infinite lattice. Our computations are performed using a 
numerical scheme based on the anti-continuous limit and Newton method 3 . 

Test and range of validity. First, we have computed numerically 
SAB (i.e. breathers whose amplitudes go to zero when Wb — > 2 + ) in the 
case B > 0. We have obtained breathers with symmetries u n (t) = u_„(t) 
(Page mode) and u n (t) = U- n -i(t + Tf,/2) (Sievers-Takeno mode), where 
Tb = 2n/ujb is the breather period. The force y n — V'{u n ) is the variable 
used in reference 2 . In Fig.l (left) it is shown that the maxima of the force 
are of order /i 1 / 2 when /i = Wb — 2 — * + , as predicted by the theory, up to 
relatively large values. Thus if B > breathers exist for any small value of 
energy in our FPU system (1). 
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Another property of these SAB is that their width diverges when Wb — > 
2+ . More precisely the theory predicts that their spatial extend is of order 
jjT 1 / 2 , which is in accordance with our numerical observations. 

Other numerical observations. For B > 0, we have numerically 
continued the SAB as ujb goes away from the phonon band. We have found 
that the maximum amplitudes of the oscillations are also approximately 
linear functions of /i 1 / 2 . This is expected for small fi, since u n = yn + 0(y 2 ), 
but it occurs surprisingly far from the phonon band, at least until values 
of fM w 1 (see fig.l, left). We have also checked that the Page mode fits 
very well to the NLS solitonu„(t) = a^/JI(—l) n cos (ubt) [cosh ([3 ^fji n)] _1 , 
even far from the top of the phonon band. 




Figure 1. Left: Force (squares) and amplitudes (circles) maxima versus /j 1 / 2 . The 
cubic coefficient in V is a = —0.1 (B = 2.64). Right: Comparison between a SAB (full 
circles) for a = —0.1 and a LAB (blank squares) for a = —1/3 (B = —1) having the same 
frequency Wt = 2.01. The dashed line represents the linear phonon with frequency 2. 



For B < and V strictly convex (-^= < \a\ < -^), breathers exist 
near the top of the phonon band but they are large amplitude breathers 1 
(LAB), i.e. their amplitudes do not go to zero when wi, — > 2+. As a 
consequence there is an energy gap for breathers creation in these FPU 
systems. In figure 1 (right) we compare a SAB and a LAB having the same 
frequency Wb = 2.01. We have found LAB with the same symmetries as 
SAB (Page and Sievers-Takeno modes) . The Page mode fits very well to an 
exponential profile having the form u n (t) = a{u>b) {— 1)™ cos (u>bt) |cr(wf,)|' ra ' 
where cr(w 6 ) = l-(^)/2 + (w fe /2)(a;^-4) 1 / 2 g (-1,0). As (1) is formulated 
as a mapping in a loop space 2 and uJb > 2, the linearized operator has a 
purely hyperbolic spectrum and the constant a(iVb) is the closest eigenvalue 
to —1 (with <r(2) = —1). Consequently, for uob ~ 2 one can ask if the iterated 
map admits a global center manifold containing these LAB. 
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